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1.   INTRODUCTION 

In  spite  of  its  familiarity  to  mankind,  boiling  liquids 
have  not  been  given  much  consideration  until  very  recently. 
The  tremendous  study  of  boiling  stems  from  the  development  of 
high  heat  flux  equipment  such  as  nuclear  reactors  and  rocket 
engines.  The  heat  flux  carried  by  boiling  is  usually  the 
order  of  106  BTU  per  square  foot  per  hour,  which  is  far  be- 
yond what  a  non-boiling  liquid  can  achieve. 

1-1.  Review  of  Previous  Developments 

The  general  subject  of  boiling  heat  transfer  is  dis- 
cussed in  detail  in  references  (l)-(6)*.  The  systematic 
study  of  boiling  began  in  1934  when  Nukiyama  (7)  published 
his  first  boiling  curve.  Since  then  rather  wide  attention 
has  been  received.  Prior  to  1950,  all  the  studies  were 
limited  to  pure  empirical  correlation.  Some  equations  had 
been  published  (8-12),  but  no  one  of  them  was  successful 
enough  to  warrant  widespread  adoption  (4). 

In  1951,  a  semi-theoretical  method  was  adopted,  which 
resulted  in  the  publication  of  the  well-known  Rohsenow's 
equation  (13)  and  the  equation  of  Forster  and  Zuber  (14). 
Both  assumed  the  major  portion  of  the  heat  is  transferred 
from  the  solid  to  the  liquid  bulk;  the  contribution  of  bubble 
latent  heat  was  considered  negligible.  This  was  based  mainly 


"Numbers  in  parenthesis  refer  to  references  listed  at  the 
end  of  the  report. 


on  the  experiments  (15-16)  observed  from  subcooled  liquids 
where  the  visible  latent  heat  carried  by  bubbles  contributed 
only  1%   or  2%   of  the  total  heat  flux.  Heat  transfer  is  con- 
sidered as  turbulent  forced  convection  stirred  by  bubble 
growth  and  detaching  velocities.  Thus,  most  of  the  proposed 
correlations  (13)  (14)  (17)  (IS)  were  of  the  form 

N   -  Const  (NRe)m  (Np  )n 
Nu 

where  the  constant  and  the  exponents  m  and  n  are  determined 
from  experiments.  At  a  particular  pressure  and  for  a  given 
surface-liquid  combination,  this  equation  reduces  approxi- 
mately to 

q/A  =  Const  (Atsat)a 

But  this  form  of  correlation  fails  to  take  care  of  the 
nucleate  characteristics  of  the  heating  surface.  Then 
Yamagata  and  Nishikawa  (19)  proposed  a  revised  form.  That 

is: 

q/A  =  Const  Nb  (Atsat)C 

1 

The  exponents  determined  by  Yamagata  and  Nishikawa  were  b  =  4, 

3  —  ■£• 

c  =  2;  by  Zuber  (20)  were  b  =  3,  c  =  3;  and  by  Tien  (21)  were 

b  =  2,  c  =  1.  A  more  detailed  discussion  of  this  development 
was  treated  by  Sato  (22)  and  Zuber  (23). 

Up  to  the  end  of  the  1950Ts,  it  was  concluded  that  latent 
heat  transport  played  only  a  minor  role  in  nucleate  boiling. 
But  this  conclusion  has  had  considerable  modification  since 


Moor  and  Mesler  (24)  have  successfully  measured  the  unusual 
cooling  effect  of  the  bubble  base.  At  about  the  same  time, 
Bankoff  (25)  carried  out  another  experiment  which  showed  that 
high  heat  transfer  rates  existed  at  the  boundary  between  the 
steam  bubble  and  the  turbulent  subcooled  liquid  stream.  A 
critical  survey  of  this  advancement  was  also  done  by  Bankoff 
(26)  in  1962.   Later  a  further  report  of  the  bubble  cooling 
effect  was  published  by  Rogers  and  Mesler  (27).  The  most 
recent  report  about  latent  heat  transport  is  due  to  the  work 
of  Rallis  and  Jawurek  (28).  They  showed  that  latent  heat 
transport  (q/A)LH  is  significant  at  all  stages,  the  ratio 
(q/A)LH/(q/A)T0T  increases  steadily  with  increasing  heat 
flux  and  appears  to  tend  to  unity  as  the  total  heat  flux 
tends  toward  burnout. 

Now  the  importance  of  latent  heat  transport  seems  rather 
clear.  In  general,  the  total  observed  heat  flux  is  contri- 
buted by  microlayer  convection  together  with  latent  heat 
transport  and  the  latter  becomes  dominant  as  heat  flux  ap- 
proaches the  burnout  point.  But  to  the  writer's  knowledge, 
no  correlation  based  on  the  latent  heat  transport  has  yet 
been  attempted. 

1-2.   Purpose  and  Outline  of  the  Report 

Although  it  seems  clear  that  latent  heat  transport  in 
nucleate  boiling  is  significant,  further  verification  is 
required.  The  purpose  of  this  report  is  to  substantiate 


the  heat  flux  contributed  by  latent  heat  transport  and  to  give 
a  theoretical  analysis  of  this  mechanism. 

In  Section  2,  some  fundamental  features  of  boiling  heat 
transfer  are  summarized  which  will  be  helpful  in  the  develop- 
ment of  the  main  part  of  this  report. 

In  Section  3,  typical  data  of  saturated  boiling  heat 
transfer  are  analyzed  first  to  show  the  numerical  fraction  of 
the  visible  latent  heat  contributed  in  saturated  liquid  boil- 
ing. Then  the  path  of  the  heat  flow  passing  through  the  bub- 
ble is  studied.  A  triple  interface  evaporation  mechanism  is 
thus  proposed.  Finally,  the  peak  flux  and  the  affecting 
factors  are  discussed. 

1-3.  Significance  of  the  Results 

A  brief  calculation  from  the  data  of  Westwater  and 

Santangelo  (29)  showed  that  the  visible  bubble  latent  heat 

carried  by  a  detaching  bubble  is  about  50%   of  the  total  heat 

flux  observed.  Careful  examination  reveals  that  the  heat 

flux  passes  through  the  bubble  via  a  very  slim  area  near  the 

solid-liquid- vapor  triple  interface  (or  line).   The  triple 

interface  evaporation  mechanism  proposed  is  based  on  this 

investigation.  Following  this  mechanism,  bubble  growth  rate 

is  given  as 

B  i 
r  =  Const  (G  tan  2)2 


which  agrees  with  the  recent  experiment  observed  by  Johnson 


and  others  (30).  An  analytically  derived  correlation  equa- 

tion  has  the  form  L 

q/A  -  Const  Nd0  — y—  exp  (-  Rf~~")Atsat 
RT|at         sat 

where  the  constant,  according  to  kinetic  theory  of  gas  and 
liquid,  should  depend  on  the  type  of  liquid.  An  explanation 
of  the  affecting  factors  using  this  equation  has  worked  very 

well. 

Since  the  correlation  made  in  this  report  is  due  to  the 
part  of  heat  flux  contributed  by  latent  heat  transport  only, 
it  can  be  considered  correct  only  at  the  high  heat  flux 
range  where  convective  heat  transfer  is  negligible.  How- 
ever, the  equation  for  bubble  growth  rate  is  always  correct. 


2.   SOME  FUNDAMENTAL  FEATURES  OF  NUCLEATE 
BOILING  HEAT  TRANSFER 

2-1.  Regimes  of  Boiling  Heat  Transfer 

The  systematic  study  of  boiling  heat  transfer  began  with 
the  discovery  of  the  unstable  region  by  Nukiyama  (7)  in  1934. 
While  boiling  a  pool  of  water  with  an  electrically  heated 
wire,  Nukiyama  found  several  regimes  existing  as  the  temper- 
ature driving  force  increased  gradually.  Fig.  1  shows  these 
typical  regimes  in  pool  boiling.   In  region  AB,  though  the 
wall  temperature  is  higher  than  the  saturation  temperature, 
no  boiling  occurs  at  the  heating  surface.  Water  evaporates 


u. 
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A t sat  =tw-tsat  ,  °F 

Fig.  1,      Regimes    oi    pool     boiling 


only  from  the  free  surface  of  the  liquid.   The  condition  is 
exactly  the  same  as  natural  convection;  heat  flux  q/A  is 
proportional  to  At*A.   In  region  BC,  bubbles  form  at  the 
active  sites  on  the  heating  surface  and  rise  through  the 
pool;  the  heat  flux  increases  very  rapidly  as  the  temper- 
ature increases;  and  q/A  is  roughly  proportional  to  At*, 

*         o   *«  A   This  region  is  called  nucleate 
where  n  ranges  from  2  to  o.   mis  region 

boiling.  At  point  C,  the  heat  flux  goes  through  a  maximum, 
after  which  q/A  decreases  as  At  increases.  This  region,  CD, 
is  called  the  unstable  or  transition  region,  and  C  is  the 
peak  point  or  burnout  point.  The  drop  of  heat  flux  in  range 
CD  is  because  most  of  the  heating  surface  is  covered  with 
vapor  film  which  prevents  the  conduction  of  heat  directly 
to  the  liquid.  At  a  point  near  D,  the  heating  surface  is 
already  completely  covered  by  a  vapor  film;  hence  region  DEF 
is  known  as  film  boiling. 

Among  all  of  these  regions,  the  nucleate  boiling  is  most 
significant  here  because  it  transfers  a  large  amount  of  heat 
in  a  moderate  temperature  difference.   The  problem  in  nucleate 
boiling  is  so  complicated  that  even  after  30  years  of  study, 
a  correct  mechanism  from  which  some  of  the  heat  transfer 
problems  can  be  predicted  still  has  not  been  devised  (31). 

2-2.   General  Consideration  of  Bubble  Dynamics 

The  unusual  heat  flux  in  nucleate  boiling  is  generally 
believed  to  be  the  result  of  the  ebullition  of  bubbles  from 


the  heating  surface.  Consequently,  to  discuss  the  phenomena 
of  nucleate  boiling  heat  transfer,  it  is  convenient  to  begin 
with  a  survey  of  bubble  dynamics. 

(a)  Active  Sites 

Liquids  have  a  tendency  to  evaporate  when  their  intrinsic 
vapor  pressure  is  greater  than  the  surrounding  pressure.  But 
this  is  true  only  when  the  liquid-vapor  interface  is  flat. 
On  a  curved  interface,  as  in  the  case  of  a  bubble,  Gibbs  (32) 
showed  that  the  equilibrium  pressures  are 

(Py  -  P)  =  R-l    R2  {1} 

which  reduces  to 

—  (9) 

(P  -  P)  =   R  UJ 

v 

for  a  spherical  bubble  as  proved  by  Kelvin  in  1870. 

From  Eq.  (1),  it  is  apparent  that  forming  a  bubble  from 
nothing  (zero  radius)  would  require  an  infinitely  large  vapor 
pressure.  It  also  means  that  only  an  infinitely  superheated 
liquid  can  ebulliate  a  bubble  from  a  perfectly  flat  surface 
with  completely  distilled  liquid.  Hsu  (33)  has  carried  out 
the  boiling  of  pure,  degassed  water  on  a  thin  layer  of  clean 
mercury  and  showed  that  no  boiling  occurred  except  an  irreg- 
ular explosive-like  formation  of  a  huge  bubble  due  to  cavi- 
tation.  This  perfect  condition  seldom  exists  in  actual  cases. 
In  ordinary  equipment,  because  of  the  nucleation  sites  on  the 


heating  surface,  bubbles  are  generated  at  a  moderate  super- 

heat  of  30F  or  less. 

Viewing  the  active  sites  directly  with  a  microscope, 
H.  B.  Clark  and  others  (34)  have  observed  that  these  sites 
are  sharp,  deep  pits  or  scratches  with  diameters  of  about 
0.0003-0.003  inch.  It  is  also  generally  believed  that 
these  pits  or  scratches  must  contain  trapped  gas  if  they 
are  to  be  capable  of  causing  bubble  generation  (35). 

(A)  O)       (C)\  (D) 

Rg.2.  Typical  shapes  of  cavity 

Consequently,  shapes  (C)  and  (D)  in  Fig.  2  are  excellent 
active  sites. 

(b)  Bubble  Stability 

From  Kelvin's  equation  (2),  it  can  easily  be  seen  that 
under  a  constant  pressure  process,  the  process  is  unstable 
(36).  Once  equilibrium  is  established,  an  accidental  in- 
crease of  radius  will  reduce  the  required  equilibrium  pres- 
sure difference.   If  the  pressure  difference  remains  un- 
changed, the  surplus  vapor  phase  pressure  will  cause  a  fur- 
ther increase  in  bubble  radius  together  with  a  further  reduc- 
tion of  the  required  balancing  pressure.  This  repeated  action 


10 


certainly  will  make  the  growth  continue  without  end.  If  a 
bubble  contains  only  inert  gas,  an  increase  of  radius  will  be 
followed  by  a  decrease  of  pressure  and  the  bubble  will  cease 
to  grow  to  an  appropriate  size.  For  a  vapor  bubble,  the  pres- 
sure of  the  vapor  inside  is  maintained  by  continuous  evapora- 
tion from  the  surrounding  liquid.   Therefore,  a  growing  bubble 
will  keep  growing  as  long  as  the  liquid  temperature  is  main- 
tained. Similarly,  a  collapsing  bubble  will  continue  collap- 
sing. For  later  reference,  two  typical  cases  will  be  discus- 
sed  below. 

(i)  Bubbles  attached  to  the  heating  surface. 

As  has  been  discussed,  once  a  nucleus  is  formed,  the 
vapor  bubble  becomes  unstable.  So,  if  the  bulk  tempera- 
ture is  high,  as  shown  in  Fig.  3  (a),  the  bubble  will 
keep  growing  until  a  certain  size  is  reached  at  which 
time  it  detaches  from  the  surface.   If  the  bulk  temper- 
ature is  not  high  enough,  as  in  the  highly  subcooled 
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Fig.  3.   Subtle  CjroM  on  the  ujali 
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case,  the  saturation  temperature  thickness  is  less  than 
the  detaching  radius,  as  shown  in  Fig.  3  (b).   The  bub- 
ble will  first  grow  to  a  size  at  which  the  vapor  condensed 
in  the  top  is  more  than  that  evaporated  from  the  bottom. 
When  this  stage  is  reached,  the  bubble  collapses.  Due 
to  the  inertia  of  the  growth,  the  bubble  will  oscillage 
in  the  range  between  R±   and  R£  as  can  be  seen  in  a  kit- 
chen kettle.  Generally  it  can  oscillate  several  times 
without  detaching  from  the  surface. 

(ii)  Bubbles  floating  in  an  infinite  medium. 

When  a  bubble  is  floating  in  bulk  liquid,  it  will 
continue  to  grow  and  absorb  some  heat  from  the  liquid, 
if  the  liquid  is  in  a  superheated  condition.  But  in 
ordinary  cases,  the  bulk  liquid  seldom  is  superheated. 
If  it  is,  the  degree  of  superheat  is  generally  limited 
to  IF  or  2F.  Hence  it  is  reasonable  to  say  that  a  bub- 
ble collapses  more  or  less  in  the  bulk  liquid. 

(c)  Bubble  Growth  Rate 

In  order  to  study  the  stirring  effect  of  bubbles,  some 

considerations  about  bubble  growth  rate  in  the  liquid  are 

necessary.   The  first  consideration  is  the  modified  Rayleigh 
equation 

R  df§  +  3  (|R,2  +  2|  .  Vll  (3) 
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where  a  term  of  surface  tension  was  added  to  his  original 
equation  (37).   Before  this  equation  can  be  solved,  some 
supplementary  relations  must  be  established.  First,  since 
most  of  the  experiments  are  based  on  temperature,  the  re- 
lation between  (Pv  -  P)  and  (tw  -  tsat)  must  be  clarified. 
In  a  superheated  liquid,  F,  Romie  (38)  has  reported  that  the 
pressure  variation  in  the  Rayleigh  equation  should  be  ob- 
tained from  the  Clausius-Clapeyron  equation 

L  L 

(Py  -  £.)  =AP  =  T(VX  -  V2)4T  =  TVfg   (tw  -  tsat) 

Second,  the  heat  of  vaporization  is  due  to  conduction  from 
liquid  bulk  to  bubble  surface.  So  there  must  be  a  relation 
between  rate  of  heat  conduction  and  bubble  radius.  This  can 
be  easily  shown  as 

dV  ^> 

L?vdS=     MAt 

k         "5      M  2   ^ 

Or,   for  a  spherical  bubble,   V  -  3  W,   d9  -  4TTR     d0,   hence 

dR     =  hto  (5) 

d0         L^ 

By  using  these  relations,  Plesset  and  Zwick  (39)  obtained  a 

relation, 

r  «/J 2KAt_ei  (6) 

which  predicts  that  the  radius  increases  as  the  square  root 
of  time.  Apparently,  Eq.  (6),  which  was  developed  for  bubble 
growth  in  an  infinite  medium,  cannot  be  regarded  as  correct 
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for  bubble  growth  on  a  wall.  Though  many  authors  have  pro- 
posed a  variety  of  approaches,  they  are  too  complicated  to 
be  considered  in  this  fundamental  treatment.  A  much  simpler 
treatment  based  on  a  new  mechanism  will  be  discussed  in  Sec- 
tion 3  of  this  report. 

2-3.  Factors  Affecting  Nucleate  Boiling 
Heat  Transfer 

(a)  Nature  of  the  Surface 

Jakob  and  Fritz  (1)  found  that  a  grooved  copper  surface 
adsorbed  air  and  initially  gave  a  much  higher  coefficient  in 
the  range  of  moderate  heat  flux.  However,  with  continued 
boiling  the  coefficient  decreased  and  very  closely  approached 
those  for  smooth  chromium-plated  surfaces.  In  both  cases 
the  area  was  taken  as  the  projected  area  of  the  plate,  re- 
gardless of  the  fact  that  the  grooved  surface  had  an  actual 
area  1.8  times  that  of  the  smooth  plate.  Deutsch  and  Rhode 
(40)  boiled  distilled  water  at  atmospheric  pressure  with 
high  heat  flux  and  found  that,  for  a  given  At,  the  coef- 
ficient U  based  on  the  projected  area  was  not  increased  by 
roughening  the  surface  and  that  U  based  on  the  total  surface 
was  less  than  that  for  the  smooth  tube. 

From  the  above  experiments,  it  can  be  concluded  that  an 
artificially-roughened  surface  will  shift  the  moderate  heat 
flux'  part  of  the  boiling  curve  based  on  the  projected  area  to 
the  low  At  end  while  maintaining  the  high  heat  flux  part  prac 
tically  unchanged  as  in  the  smooth  surface. 
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For  surfaces  of  different  materials,  the  data  of  Bonilla 
and  Perry  (12)  (Fig.  4)  show  that  for  boiling  ethanol  at  at- 
mospheric pressure  on  a  horizontal  flat  plate,  a  higher 
thermal  conductivity  gives  higher  heat  flux  at  a  given  At, 
while  the  maximum  heat  flux  remains  practically  the  same  for 
all  four  surfaces.  But  with  water  or  methanol  boiling  at 
atmospheric  pressure  in  a  small  submerged  tube  evaporator, 
Cooper  and  others  (40)  found  that  U  was  larger  with  iron  than 
copper  tubes,  indicating  that  the  increase  in  the  number  of 
vaporization  nuclei  had  over-compensated  for  the  decrease  in 
thermal  conductivity.  Thus  it  can  be  concluded  that,  for 
different  surfaces,  the  higher  the  thermal  conductivity  and 
the  more  the  evaporation  nuclei,  the  higher  the  resulting 
heat  flux.  An  aged  surface  generally  decreases  both  thermal 
conductivity  and  number  of  evaporation  nuclei. 

(b)  Effect  of  increased  velocity 

In  the  range  of  low  At,  the  use  of  forced  convection  in 
a  boiling  system  results  in  an  increase  in  the  heat  flux  for 
a  given  At.  But  in  the  region  of  strong  nucleate  boiling,  the 
influence  of  velocity  is  small,  as  shown  in  Fig.  5.  Line  A 
shows  data  of  Beecher  (40)  for  water  at  212F  and  1  atm  flowing 
at  3  ft./sec.  normal  to  an  electrically  heated  0.050-inch   . 
diameter  stainless  steel  tube;  Line  B  represents  data  for 
water  boiling  on  a  0.046  inch  platinum  wire  in  an  unstirred 
pool  at  1  atm.  It  is  apparent  from  this  figure  that  in  a 
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Fig. 4.    Effect    of   heating  surface  02) 


strong  nucleate  boiling  range  the  heat  flux  is  essentially 
unaffected  by  mechanical  agitation. 
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(c)  Effect  of  subcooling 

In  an  ordinary  heat  transfer  problem,  the  heat  flux  is 

generally  proportional  to  the  temperature  driving  force 

(t  -  t  ),  but  things  are  always  unusual  whenever  boiling 

x  w    "° 

heat  transfer  is  concerned.  Figure  6  shows  a  logarithmic 
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graph  of  the  heat  flux  q/A  plotted  as  ordinate  vs.  the  total 
Z\t  from  heater  to  the  degassed  water.   In  the  nonboiling 
region,  the  results  for  each  agree  with  those  expected  from 
conventional  equations  for  forced  convection  without  change 
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in  phase.   In  the  local  boiling  region,  the  curves  are  steep 
and  are  displaced  horizontally  for  each  value  of  subcooling 
by  values  of  At  corresponding  closely  to  the  differences  in 

subcooling. 

When  the  same  data  for  surface  boiling  are  plotted,  in  a 

different  fashion,  as  shown  in  Fig.  7,  with  Atgat  "  *w  -  tsat 
replacing  At  -  tw  -  t^  as  abscissa,  the  results  are  insensi- 
tive to  water  temperature  and  velocity.  The  slope  of  the  curve 
in  Fig.  7  is  similar  to  that  for  the  boiling  of  a  pool  of 
saturated  liquid. 

(d)  Effect  of  pressure 

For  a  large  number  of  liquids  boiling  in  pools  in  the 
nucleate  region  at  pressures  of  1  atm.  and  less,  the  data 
show  that  a  decrease  in  saturation  pressure  gives  a  lower 
heat  flux  for  a  given  At  (40).  For  pressures  higher  than 
atmospheric  pressure,  the  data  also  show  that  an  increase 
in  saturation  pressure  gave  a  higher  heat  flux  for  a  given 

At,  but  the  peak  flux  is  a  curved  line  with  a  maximum  mid- 
way to  the  critical  state.   Figure  8  shows  the  curve  plotted 
from  the  data  of  Addom  (3).  The  heat  transfer  coefficient  at 
2465  psia  was  100  times  greater  than  the  value  for  the  same 

At  at  1  atm. 

2-4.   Some  Proposed  Mechanisms 

Several  mechanisms  of  nucleate  boiling  heat  transfer 
have  been  proposed  during  the  past  three  decades.  They  were 
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all  suggested  to  correlate  the  data  of  the  extremely  high 
heat  flux  of  nucleate  boiling.  Unfortunately  no  one  of  them 
had  ever  successfully  explained  the  experimental  facts  stated 
in  the  preceding  section. 
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Fig.  6.     If  fed  of    pressure   (3) 


(a)  Microconvection  in  the  sublayer 


This  mechanism  is  most  widely  accepted  at  the  present 
time.   The  heat  path  is  assumed  to  lead  from  the  heating  sur- 
face  to  the  liquid  between  the  bubbles.   Observed  high  heat 
flux  is  considered  due  to  the  stirring  effect  of  the  bubbles, 
As  shown  in  Fig.  9,  contrary  to  the  convective  velocities 
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Fig.  9.      Stirring  effect  of  bubbles 


that  change  the  sublayer  pattern  from  outer  layers,  the  bub- 
ble growth  velocities  stir  the  liquid  from  inside  the  sub- 
layer.  Consequently,  these  radial  velocities  of  the  order 
of  10  or  20  fps  are  estimated  so  large  that  they  determine 
the  temperature  distribution  in  the  sublayer  near  the  heat- 

ing  surface. 

This  mechanism,  of  course,  provides  a  high  heat  flux  in 
nucleated  boiling  if  everything-  happens  to  be  as  described 
above.   However,  it  is  also  very  important  to  observe  that 
if  heat  flows  from  the  heating  surface  through  the  thermal 
sublayer  to  the  bulk  liquid,  it  depends  strongly  on  the 
temperature  difference  (tw  -  t„  )  which  is  the  driving  poten- 
tial of  the  heat  flux.  .  But  as  already  discussed  in  Section 
2-3,  the  heat  flux  is  essentially  independent  of  the  degree 
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of  subcooling.  This  proposed  mechanism  is  not  as  good  as 

originally  expected. 

(b)  Bubbles  act  in  the  manner  of  surface  roughness 

That  bubbles  act  in  the  manner  of  surface  roughness  was 
suggested  by  H.  S.  Tsien  (3D-  For  boiling  heat  transfer 
with  forced  convection,  the  bubbles  on  the  heating  surface 
increase  the  turbulent  exchange  of  liquid  between  the  heating 
surface  and  the  moving  bulk  liquid.   The  effect  is  similar  to 
that  produced  by  plate  roughness.   But  just  as  for  the  mech- 
anism (a),' the  heat  flux  would  again  have  to  depend  directly 
on  the  temperature  driving  force  tw  -  t„  .  Furthermore,  the 
turbulent  exchange  has  nothing  to  do  with  the  case  of  pool 
boiling,  which,  it  is  believed,  should  have  the  same  basic 
mechanism  as  in  forced  convection  boiling. 

(c)  Latent  heat  transported  by  bubbles 

Latent  heat  transport  means  that  while  a  bubble  grows  it 
absorbs  the  latent  heat  of  vaporization  which  is  then  returned 
to  the  bulk  liquid  where  the  bubble  collapses.  By  calculating 
the  total  bubbles  formed  and  the  total  latent  heat  contained 
in  these  bubbles,  Rohsenow  and  Clark  (15),  also  Gunther  and 
Kreith  (16)  showed  that  the  latent  heat  carried  by  the  bubble 
contributed  only  1%   or  2%   of  the  total  heat  flux.   It  is  in- 
teresting to  note  that  both  cases  studied  were  under  conditions 
of  high  subcooling  of  150F  and,  as  pointed  out  by  Bankoff  (26) 
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and  Snyder  (41),  this  heat  flux  carried  by  latent  heat  could 
be  increased  by  assuming  additional  heat  flux  through  the 
bubbles  by  mass  transfer.  That  is,  heat  flows  into  the 
individual  bubble  through  the  superheated  base  near  the  heat- 
ing plate.  This  heat  is  absorbed  as  heat  of  vaporization  at 
the  vapor-liquid  boundary,  and  is  then  carried  as  steam  to 
the  top  of  the  bubble  where  the  vapor  condenses  and  gives  off 
latent  heat  to  the  subcooled  liquid  bulk.  A  more  detailed 
discussion  will  be  treated  in  Section  3  of  this  report. 

(d)  Vapor-liquid  exchange  action 

This  vapor-liquid  exchange  action  was  a  mechanism  sug- 
gested by  Forster  and  Greif  (U).  The  main  idea  is  that  when 
a  bubble  is  formed  and  detaches  from  the  heating  surface  or 
collapses  above  the  heating  surface,  it  pushes  an  amount  of 
hot  liquid  having  the  same  volume  as  the  bubble  into  the 
liquid  bulk,  as  shown  in  Fig.  10.  This  vapor-liquid  exchange 
is  assumed  to  contribute  the  most  heat  flux  in  nucleate 
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Fig.  10.     Vapor  liquid  exchange  action 
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boiling  heat  transfer.   Indeed,  by  using  the  following 
equation,  Forster  and  Greif  (41)  were  able  to  show  that  the 
heat  transferred  by  this  vapor-liquid  exchange  action  would 
contribute  about  320$  of  the  total  observed  heat  flux. 


AH  =  C  9   V  At 


(7) 


mean 


But  this  only  means  that  vapor-liquid  exchange  action  could 
be  a  main  contributer  to  nucleate  boiling  heat  transfer  only 
if  the  temperature  difference  Atmean  used  in  Eq.  (7)  is  cor- 
rect or  nearly  correct.  From  Fig.  11,  it  is  easily  seen  that 
when  a  bubble  detaches  from  the  surface  or  collapses  over  the 
surface,  the  liquid  pushed  into  the  bulk  is  actually  the  re- 
latively low  temperature  part.  Therefore,  the  temperature 
difference  to  be  used  in  Eq.  (7)  should  be  much  less  than 

At    ,  and  the  contribution  by  this  exchange  should  be 

mean' 
much  less  than  that  calculated  by  Eq.  (7). 


Liquid  exchotoged 


Fig.  II.      Temperature  of  1he  liquid  exchanged 
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(e)  Microlayer  evaporation 

The  most  recently  proposed  mechanism  is  microlayer 
evaporation,  as  suggested  by  Moore  and  Mesler  (24).  By 
using  a  sensitive,  plated  thermocouple  directly  on  the  bub- 
ble  site,  Moore  and  Mesler  observed  an  unusual  temperature 
drop  on  the  surface  while  the  bubble  was  growing.  From  this 
they  concluded  that  a  microlayer  exists  inside  the  bubble, 
and  the  evaporation  of  this  microlayer  will  contribute  a 
large  amount  of  the  heat  flux  observed.   But,  as  pointed  out 
by  Lyon,  Fourst  and  Katz  (42),  if  wetting  does  not  occur, 
then  the  formation  of  a  microlayer  would  be  unlikely.  Thus, 
nucleate  boiling  would  not  be  expected  if  it  is  dependent 
upon  microlayer  vaporization. 


Fig.  12.     Microlayer     evaporation 
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3.  TRIPLE  INTERFACE  EVAPORATION 

Investigators  constantly  point  out  that  the  latent  heat 
transport  should  have  a  major  role  in  nucleate  boiling  heat 
transfer.   However,  neither  an  experimental  method  nor  an 
analytical  approach  has  ever  been  tried  to  give  a  more  con- 
crete evaluation  of  this  process.  To  make  things  clear,  an 
analytical  study,  though  brief  in  itself,  is  necessary. 

3-1.  Contribution  of  Bubble  Latent  Heat 

Gunter  and  Kreith  (16)  and  Rohsenow  (15)  have  observed 
that,  when  boiling  a  highly  subcooled  liquid,  the  latent  heat 
contained  in  the  bubbles  has  an  order  of  magnitude  of  only 
1%   or  2%   of  the  total  observed  heat  flux.  This  does  not 
necessarily  mean  that  the  heat  transferred  to  the  bubble  is 
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Fig.  13.    Heat  flux  through  latent   heat  transport 
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insignificant.  It  can  be  seen  in  Fig.  13  that  when  boiling 
a  highly  subcooled  liquid,  vapor  bubbles  are  always  larger 
than  the  thickness  of  thermal  boundary  layer  (31).  This  ex- 
plains why  bubble  diameters  are  smaller  in  a  subcooled  con- 
dition than  in  a  saturated  liquid.   However,  due  to  the  high 
thermal  diffusivity  in  liquid  bulk,  the  thermal  layer  on  top 
of  the  bubble  can  be  expected  to  be  very  small.  Across  this 
extremely  thin  thermal  layer,  the  temperature  driving  force 

is  Wap  "  hiq  =  teat  "  *- "  degree  °f  subC0°ling*  C°RSe~ 
quently,  the  latent  heat  contained  in  the  bubble  itself  may 

be  of  a  small  order  of  magnitude.  The  heat  flux  transferred 
to  the  bulk  liquid  through  the  top  of  the  bubble  in  a  highly 
subcooled  condition  is  undoubtedly  of  very  large  magnitude. 
To  check  the  possibility,  Fig.  13  shows  that  if  heat  flux 
passing  through  B  is  very  large,  the  heat  flux  through  A 
-into  the  bubble  should  also  be  very  large.  Because  the  tem- 
perature difference  is  the  same  (tw  -  t^)  for  both  subcooled 
and  saturated  cases,  the  heat  flux  passing  through  A  should 
remain  practically  the  same  for  both  cases.   (This  can  be 
considered  as  a  reasonable  assumption  for  the  time  being  and 
will  be  proved  later).  Furthermore,  in  the  saturated  case 

t         „  t     ■  t    -t    =0,  the  heat  flux  through  B 
^vap    bulk    sat    sat 

should  be  zero.  All  the  heat  flux  through  A  must  equal  the 
latent  heat  of  vaporization  of  the  vapor  contained  in  the 
bubble.   This  indicates  that  if  one  can  find  a  large  amount  of 
heat  carried  by  the  vapor  contained  in  the  bubble  when  boiling 
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saturated  liquid,  it  can  be  concluded  that  the  latent  heat 
transport  at  the  top  of  bubble  is  significant  when  subcooled 
liquid  is  boiled.  This  is  the  actual  case  if  a  brief  calcu- 
lation is  made  from  the  data  of  Westwater  and  Santangelo  (29). 

(A)  Data  of  Westwater  and  Santangelo: 
.  Liquid:  methanol  at  1  atm 

.  Observed  overall  coefficient:   U  =  1350  BTU/hr  ft2F 

.  Overall  temperature  difference:  At  -  70F 

.  Average  bubble  size:   d  =  0.17  inch  diameter 

.  Bubbling  rate:      f  =  17  bubbles/sec. 

.  Spacing  of  nuclei  sites  =  0.103  inch 

(B)  A  calculation  based  on  one  square  foot: 

.  Total  nuclei  site  =  ^-^§j  x  ^-^|j  =  H400 

.  Total  number  of  bubbles  leaving  heating  surface 

per  hour  =  14400  x  17  x  60  x  60  =  3.8  x  10 
.  Volume  of  one  bubble  «  3>14  *  °'^   =  0.00257  in3 

-  1.485  x  10"  ft3 
Specific  volume  of  CH30H  at  saturation  (43) 

-  13.05  ft3/lb  _6 

.  Mass  of  one  bubble  =  1'k%%°       "   1.138XHT7  lbs 

.  Latent  heat  of  CH  OH  =  482  BTU/lb  at  1  atm. 

.  Total  heat  carried  by  bubbles 

-  482  x  1.138  x  10"7  x  B.B   x  10g 

=  48200  BTU/hr  ft2  =  688  BTU/hr  ft2  F 

.  Percentage  of  observed  heat  flux  carried  by  de- 
taching bubble  is   H|o-  =  ^ 
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This  is  the  precentage  contributed  by  bubble  latent  heat  in 
a  condition  still  much  below  the  peak  condition.  At  peak 
flux,  q/A  =  172,000  BTU/hr  ft2,  the  percentage  will  be  still 

more. 

The  high  heat  content  of  a  vapor  bubble  also  can  be 
easily  visualized  in  ordinary  high-duty  power  boilers.  In  a 
power  boiler  as  shown  in  Fig.  14,  all  the  heat  absorbed  from 
the  heating  surface  is  used  to  generate  saturate  steam.   If 


Fig.  14.   Schematic  diagram  of  poajer   boiler 
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the  steam  is  generated  by  evaporation  from  the  free  surface, 
then  owing  to  the  large  amount  of  heat  absorbed  in  evapora- 
tion, the  temperature  at  A  should  be  much  higher  than  the 
saturate  steam  temperature  at  B.  However,  the  temperature  at 
A  differs  only  slightly  from  tB  as  can  be  seen  from  any  boiler 
operating  data  (44).  From  this  fact,  it  is  reasonable  to  say 
that  the  steam  to  be  generated  has  already  become  steam  while 

still  inside  the  boiler  tubes.  This  certainly  means  that  the 

absorbed  heat  has  gone  entirely  to  the  bubbles. 

Although  only  the  case  of  boiling  a  saturated  liquid  is 

discussed  here,  it  will  be  proved  later  that  the  condition  is 

the  same  when  a  subcooled  liquid  is  boiled. 

3_2.  Triple  Interface  Evaporation 

The  high  heat  flux  carried  away  by  bubbles  when  boiling 
a  saturated  liquid  has  already  been  discussed.  The  next  step 
is  to  find  how  such  a  large  amount  of  heat  can  flow  into  the 
bubbles.  A  general  discussion  of  interface  heat  transfer 
follows.  From  the  kinetic  theory  of  gases  and  liquids,  it 
is  known  that  liquid  molecules  evaporate  continuously  when- 
ever the  surrounding  pressure  is  less  than  the  intrinsic   • 
vapor  pressure  of  the  liquid.   (Of  course,  only  net  evapor- 
ation is  of  interest  here.)  Also,  the  latent  heat  carried 
away  by  this  evaporation  is  so  fast  that,  when  compared  with 
the  heat  rate  conducted  from  the  liquid  bulk,  it  can  be  con- 
sidered as  infinite.  This  infinite  heat  transfer  coefficient 
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is  actually  an  assumption  adopted  by  most  authors  when  treat- 
ing the  bubble  growth  rate  in  boiling  (33)  (45).  Experi- 
mentally, Alty  (46)  also  observed  that  when  reducing  the 
vapor  phase  pressure,  a  water  droplet  temperature  of  25C 
was  able  to  form  a  thin  layer  of  ice  on  the  surface.   For  a 
numerical  concept,  it  is  known  that  the  maximum  rate  of  evapor- 
ation of  a  liquid  evaporating  to  a  vacuum  is  equal  to  the  rate 
of  the  vapor  molecules  that  would  collide  on  the  liquid  sur- 
face if  the  vapor  is  at  saturation  pressure.  More  precisely, 
it  equals  the  vapor  molecules  that  are  colliding  on  the 
liquid  surface  and  have  been  captured.  Thus,  the  maximum 
rate  of  evaporation  to  vacuum  can  be  expressed  as  (47) 

/  £M  .  ,£X 

m  =  £  f  P]  TTRT  [i5i 

where  f  is  a  factor  taking  care  of  the  molecules  that  collide 
with  but  are  not  captured  by  the  liquid.  The  value  of  f 
differs  from  liquid  to  liquid.   For  most  liquids,  however,  it 
is  very  close  to  unity  but  for  water  at  212F  it  was  shown  by 
Alty  (46)  that  f  has  a  value  of  only  0.04-  Therefore,  for 
water  at  212F,  Eq.  (8)  becomes 

f       _g_x  18  /     '       ,.    i 

m  =  0.01  P       334  x  49750  x  T    (slugs/ft    ,   hr) 

or  m  -  0.00975P  J  T    (lbs/ft2,   hr)  (9) 

For  water  at  248F  evaporating  to  a  vapor  space  of  1  atm,  the 
net  rate  of  evaporation  is: 
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nu   -  nu  -  0.00975  x  1U   (23.797  J708  -  14.7   j  m) 
=  2510  lbs/ft2,   hr 

The  latent  heat  carried  by  this  evaporated  vapor  is 

q/A  =  2510  (1150.4  -  216.45)  =  2.34  x  106  BTU/ft2  hr 

This  shows  that  for  a  temperature  difference  of  36F  at  1  at*. 
the  heat  flux  carried  by  evaporation  is  5  times  as  much  as  the 
peak  flux  observed  in  nucleate  boiling.  This  also  shows  that 
if  20  percent  of  the  area  goes  through  an  evaporation  process, 
it  will  result  in  a  heat  flux  equivalent  to  the  maximum  heat 
flux  in  nucleate  boiling. 

A  physical  picture  of  the  bubble  heat  transfer  process 
is  now  in  order.   In  treating  the  growth  rate  of  bubbles, 
Griffith  (45)  assumed  the  boundary  conditions  of  a  growing 
bubble  as  shown  in  Fig.  15  (A).   In  this  figure,  there  is  a 
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Fig.  15.      Discontinuity     in   triple    interface 
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discontinuous  point  called  the  triple  interface  (or  line)  A. 
At  this  point,  it  is  obvious  that  the  temperature  gradient 
across  the  corner  is  infinite.  Thus,  if  the  condition  is 
good  enough,  the  heat  rate  through  the  neighborhood  of  A  can 
be  infinite  in  magnitude.  To  check  the  behavior  of  this 
region,  a  simplified  stationary  two-dimensional  square  bubble 
is  assumed  as  shown  in  Fig.  15  (B).  Although  the  actual  bub- 
ble is  growing,  unsteady  and  spherical,  the  simplified  bubble 
does  provide  a  good  approximation  so  far  as  the  qualitative 
behavior  of  the  corner  is  concerned.  The  steady  temperature 
distribution  in  the  square  liquid  e-f-g-h  has  been  solved  by 
Carslaw  and  Janger  (48)  as 

t-W= it     „t0te»+,)        a 

In  order  to  calculate  the  heat  flow  to  the  vapor  phase,  the 
above  equation  must  be  differentiated  to  find  the  temperature 
gradient  at  x=0.  But  the  termwise  differentiation  is  not 
allowed  because  in  differentiating  the  above  equation  term  by 
term  results  in 

at  4Ck^^yjLCos(-^^ 

"sT"   it    0  a     a        a 

or 
^t,  _4(U-W)  §sir,hCa"yX2h+')rr  Coseck(2h-H)TT       00) 

which  does  not  converge  at  y=0.  But  a  little  rearrangement 
shows  that  Eq.  (10)  converges  uniformly  in  the  interval 
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s  ±  7  £  a,  where  5  >  0.   (See  Appendix  1. )  Thus  the  ternad.se 
differentiation  is  valid  for  this  interval. 

For  5  ±  J  ±   a,  the  heat  flow  from  liquid  phase  to  vapor 


phase  is 


Kaxl,=0  ' 


By  proper  simplification,  some  values  of  Q5.a  are  calculated  in 
Appendix  2.  The  general  behavior  of  the  triple  interface  dis- 
continuity is  shown  in  Fig.  16.  From  this  figure,  it  is  ap- 
parent that  the  heat  flow  into  the  vapor  phase  is  infinite. 
Also,  it  is  shown  that  this  heat  flux  is  mostly  contributed 
by  the  slim  area  near  the  bubble  base  that  is  in  contact  with 
the  heating  surface.  This  is,  of  course,  an  ideal  case  that 
would  not  occur  because  of  the  meaningless  infinite  heat  flux.  But 
it  is  true  that  the  heat  flux  distribution  in  an  actual  case 
should  resemble  the  distribution  of  that  shown  in  Fig.  16.    .  ' 
That  is,  the  heat  flow  into  the  bubble  occurs  mostly  over  the 
small  region  of  the  liquid-vapor-solid  triple  interface. 
Since  the  condition  near  the  triple  interface  should  be  the 
same  for  both  saturated  and  subcooled  cases  when  tw  remains 
the  same,  the  heat  flow  into  the  bubble  also  should  be  the 
same  under  the  same  tw  and  tgat. 
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Kg.  16.  General  behavior  of  triple  interface  dlscointinu'itij 


For  an  experimental  proof  of  this  extreme  heat  flux  at 
the  discontinuity,  it  is  interesting  to  take  a  look  at  Cos- 
tello  and  Redeker's  (49)  experiment.  By  using  capillary 
wicking  as  shown  in  Fig.  17-  Costello  and  Redeker  have  ob- 
served an  amazing  heat  flux  that  is  far  more  than  that  ob- 
tained by  ordinary  pool  boiling.  Although  the  authors  did 
not  state  definitely  the  cause  of  this  high  heat  flux,  a 
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Stainless  steel   tube 

0.125"  O.D. 


Fig.  IT.  Capillary  wicking  (49) 


tentative  conclusion  may  be  that  the  capillary  wicking  pro- 
cess provided  two  solid-liquid-vapor  intersection  lines  along 
the  longitudinal  direction  of  the  tube. 

Further  proof  of  this  discontinuous  behavior  is  a  very 
sharp  temperature  drop  which  should  be  observed  at  the  heating 
surface  as  the  bubble  boundary  passes  through.  This  is  due 
to  the  very  high  heat  flux  over  the  extremely  small  area  at 
the  base  of  bubble.  This  is  precisely  the  case  observed  by 
Moor,  Robers  and  Mesler  (24)  (27).  Meanwhile,  the  assumption 
of  microlayer  evaporation  seems  unnecessary. 

3-3   Bubble  Growth  Rate  and  Correlation 
of  Heat  Flux 

In  the  preceding,  the  extreme  behavior  of  the  triple 
interface  has  been  discussed.  For  a  mathematical  represen- 
tation, it  would  be  convenient  to  assume  an  effective  thick- 
ness S ,  within  which  liquid  evaporates  through  the  whole 
temperature  driving  force  tw  -  tsat,  and  the  contribution  is 
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completely  neglected  beyond  this  thickness.  Similar  to  the 
mixing  length  assumed  in  the  turbulent  flow  problems,  it  is 
hoped  that  this  effective  thickness  will  be  a  constant  or  at 
most  dependent  on  the  saturation  temperature  only.  In  this 
preliminary  treatment,  it  will  be  considered  as  constant. 
The  magnitude  of  S  for  water  is  about  ^  *0  (APP^ndix  3). 
Mow  the  instantaneous  area  of  a  bubble  can  be  written, 

(12) 


A  -  It  dc  S 


where  dc  is  the  instantaneous  contact  diameter. 

To  find  the  evaporation  rate  through  the  effective  area, 
Eq.  (9)  can  be  used.  However,  because  this  equation  involves 
two  independent  variables,  P  and  T,  another  means  of  expres- 
sion is  preferable.  From  the  kinetic  theory  of  liquids,  it 
is  known  that  the  rate  of  evaporation  is  proportional  to  the 
number  of  molecules  having  kinetic  energy  greater  than  their 
bonding  energy.  Thus,  for  a  liquid  following  the  Boltsmann 
energy  distribution,  its  evaporation  rate  to  a  vacuum  can  be 
expressed  as, 

ti°\  (13) 

n1  =  a  exp  v-^t' 
or,  since  N  6Q  =  L,  Nk  =  R,  Eq.  (13)  also  can  be  written  as 

/  J±\  (14) 

m  -  a  exp  \~  RT> 

where  "a"  is  a  proportionality  constant.   For  evaporation  to 
bubble,  the  net  rate  of  evaporation  is 
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L  L    -I 

m  -  a  [  exp  (-  RTj  -  e*P  (-  RTsat  )  J 


(15) 


In  the  practical  case,  Eq.  (15)  is  still  very  difficult  to 
apply.  WhenAtsat  is  saall  -pared  with  T^.  this  equation 
can  be  simplified  as, 


A  =  a"       «XP  t-^sat'  Atsat 
R1sat 


(16) 


which  shows  that  the  evaporation  rate  per  unit  area  is  appro*,- 
mately  proportional  to  the  first  power  of  the  temperature  dif- 
ferenoe,  but  differs  from  temperature  to  temperature.  A  word 
of  caution,  usually  Eq.  (16)  is  not  a  good  approximation  be- 
cause Atsat  may  become  very  large  in  which  case  Eq.  (15) 
should  be  used  instead. 

Consider  a  bubble  with  a  contact  angle  p  ,  shown  in 
Fig.  18.  Let  the  bubble  radius  increases  dr  during  a  time 
interval  d.9.  Since  the  evaporated  vapor  should  be  equal  in 
volume  to  the  bubble  volume  increase,  the  following  should 
result. 


Fig.  18.  Bubble  growth  on    the  wall 
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1  — fe—  t   pt     "Ut       d6  -  2ITr2  1+Cos  ft  )dr 

p-Hd  Sin/3Sa  RT2       exp(-RTsat Wsatao       "  V 


or 


dr 


SaL  exp  (-  RTsat)Atsat  tan  £ 


That  is, 


d6  2 

r  ?v  RTsat 


dr  „  b  t   £ 
d6    r  tan  2 


(17) 


where 


b  = 


*SL^sat  exp  <-ET^> 


,2 


?v  RTsat 
Integrating  Eq.  (17)  gives 


r.  j  2b  tan  f   ei  <l8> 

which  shows  that,  for  a  given  contact  angle   ,  the  bubble 
radius  is  proportional  to  the  square  root  of  time.  This  has 
been  shown  hy  many  authors  using  completely  different  mech- 
anisms, in  an  actual  case,  once  the  bubble  begins  to  grow, 
the  heating  surface  temperature  will  fall  considerably. 
Therefore,  the  actual  growth  rate  should  be  a  little  slower 
than  that  shown  in  Eq.  (18).  In  a  particular  bubble,  West- 
water  (50)  observed  that  r  is  proportional  to  $U'  .  For  a 
different  contact  angle,  it  is  found  from  Eq.  (13)  that 
although  the  relation  between  radius  and  time  remains  the 


40 


same,  the  constant  of  proportionality  increases  with  an  in- 
crease of  contact  angle.  That  is,  the  bubble  growth  rate  is 
higher  for  a  semi-spherical  bubble  than  for  a  spherical  bub- 
ble.  This  fact  has  been  observed  by  Johnson,  Jr.  and  others 
(30),  though  the  observers  have  a  different  explanation  for 
it.   For  simplicity,  the  rest  of  this  report  will  consider 
only  the  semi-spherical  bubble,  that  is  |S»V2,  or  dQ   -  d. 

Before  Eq.  (12)  can  be  applied  to  calculate  an  average 
effect  of  the  heat  flux,  the  mean  effective  bubble  diameter 
d*  must  be  found.   Since  bubble  diameter  changes  constantly, 
the  mean  effective  diameter  may  be  defined  as 

■ 

Substituting  Eq.  (18)  gives 

a*=  -^  J*2  J£B  8**19  =  -|-d0  (19) 

where  d  is  the  detaching  bubble  diameter. 

From  Eqs.  (12),  (16)  and  (19)  the  total  heat  flux  due 
to  nucleate  boiling  can  be  written  as 

q/A  -ON*/*  5rT-  exp  (-  R^T^t^  L 
sat 

2  L2  L 

or  q/A  «  3  «  TT   NdQ  ^2~  exp(-RTsat>Atsat  (20) 

sat 

Initially  the  heat  flux  seems  proportional  to  the  first  order 
0f  Atsat,  but  this  is  not  true  since  bubble  number  N  increases 
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very  rapidly  as  temperature  increases. 

Gaertner  and  Westwater  (51)  adopted  a  nickel-plating 
method  and  found  that  heat  flux  varies  approximately  with  the 
square  root  of  the  active  site  population.  Also,  the  jamming 
effect  (33)  and  the  high  rushing  velocity  of  evaporating  va- 
por in  a  high  superheat  state  causes  the  detaching  bubble 
diameter  (or  aureole  diameter)  to  decrease  with  an  increase 
of  superheat.  Figures  19  and  20  are  two  rearranged  curves 
from  Figs.  7  and  13  of  Gaertner  and  Westwater ■■•  paper. 

From  Fig.  19,  the  relation  between  the  detaching  bubble 
diameter  and  the  active  site  population  in  the  high  heat  flux 
range  is  found  as 

V  (21) 

"0  =  ?^ 

The  detaching  diameter  dQ  remains  roughly  the  same  when  heat 

flux  is  small. 

From  Fig.  20,  the  relation  between  4tsat  and  active 

site  population  is  obtained  as 

r 

~.T0.08l  (22) 

^sat  =  TN 

Then  from  Eqs.  (20),  (21)  and  (22),  the  following  equ- 
ations  result: 
(a)  At  low  heat  flux: 

2  L2  L     ..O.OSl 

q/A  -  3  aSFNd0  RT2    exp  (-RTsat)  0  N 


"sat 
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F/g. 20.     Relation  between  superheat   and    bubkie    population 
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or 


q/A  oZ     N 


l.OSl 


(23) 


which  shows  that  heat  flux  is  roughly  proportional  to  the 
active  site  population  as  observed  by  Jakob  (52). 

(b)  At  high  heat  flux: 

2       %  L2  L      0.081 

q/A  =  3  fl5nN  OT  RT2    exp  (-  RTsat)  T  N 


sat 


or 


q/A  ©C  N 


0.401 


(24) 


which  agrees  with  the  experiment  of  Gaertner  and  Westwater. 

If  the  heat  flux  is  expressed  by  the  term  Atsat  only, 
Eq.  (24)  becomes 

1 
r  -atea1,.  n- 

q/A  oC 


(   r   ' 


0.401 


That  is, 


q/A  oC  (4t   .  ) 


(25) 


sat 


for  this  particular  case. 

From  the  above  verification,  it  is  seen  that  Eq.  (20) 
agrees  very  well  with  the  experiments.   But  since  it  should 
be  expressed  in  terms  of  active  site  population  and  detaching 
bubble  diameter,  no  simple  relation  between  heat  flux  and 
At  *  can  be  deduced  without  a  knowledge  of  active  site 
variation  and  bubble  diameter  behavior.   This  fact  may  be 
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considered  as  the  cause  of  the  scattering  of  the  experimental 
results. 

3-4.  Maximum  Flux  and  Factors  Affecting 
Nucleate  Boiling  Heat  Transfer 

Equation  (20)  shows  that  at  a  given  temperature,  heat 
flux  is  proportional  to  bubble  population  N,  bubble  diameter 
d0  and  temperature  driving  force  ^t^.   Since  dQ  does  not 
change  very  much  withAtsat,  heat  flux  will  generally  increase 
as  N  increases.  But  this  is  true  only  when  each  bubble  is 
geometrically  independent.  When  N  reaches  a  value  at  which 
the  boundaries  of  the  bubble  begin  to  eclipse  each  other,  the 
effective  circumference  per  bubble  will  be  reduced  consider- 
ably. Certainly,  if  this  reduction  cannot  be  compensated  for 
by  the  increase  of  N,  the  heat  flux  will  then  decrease  in 
spite  of  the  fact  that  wall  temperature  is  kept  increasing. 
This  is  the  transition  region  defined  by  the  triple  interface 
evaporation  mechanism;  the  turning  point  is  know  as  the  burn- 
out point. 

To  find  the  peak  flux,  it  is  necessary  to  know  the  geo- 
metrical distribution  of  the  bubbles.  Artificially  the  maxi- 
mum number  of  bubbles  that  can  be  packed  without  their  boun- 
daries touching  each  other  is  the  quantity  obtained  by  Fig.  21 
(A).  This  is  the  well-known  hexagonal  packing.  The  ratio  of 
the  bubbles  obtained  by  this  hexagonal  packing  to  that  of 
ordinary  square  packing  is  a/b  -  l/sin  60°  =  1.15.  Though 
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04)      Hexagon  (B)      Square 

Fig.  21.    Bubble    packaging 


the  bubbles  as  arranged  in  Fig.  21  (A)  are  not  necessarily 
the  ones  that  give  the  maximum  heat  flux,  it  is  very  reason- 
able to  assume  that  the  maximum  heat  flux  happens  when  the 
mean  effective  diameters  are  arranged  as  shown.  From  this 
postulate,  the  following  equation  results  at  burnout  point: 


N 


t  = 


1  •  _Jk 
1.15  (d*   d*} 


1 

1.15  (4 


1 


1_ 

0 


do 


1  ,1  '  . 

2d*  =  1.0.5    (2    Jfi\) 


(26) 


Substituting  Eq.  (26)  into  Eq.  (20)  yields 

r-         _JdL  L 

(q/A)v  =  1.05  aSirVN'   RT2    exp  (-  RTsat)At^at   (27) 


max 


RT 


sat 
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or,  in  terms  of  d', 


L2  L 


<^>max  "  ^  a5,f  di   RT2at  eXp("  ^sat^sat   <28> 

Before  these  equations  can  give  the  exact  value  of  the 
heat  transferred,  a,  $  ,  d^  and  Atgat  should.be  determined 
either  experimentally  or  by  statistical  methods.  However, 
these  equations  are  already  sufficient  to  be  used  to  predict 
the  qualitative  effect  of  one  factor  upon  the  others.  Some 
of  these  effects  are  discussed  below. 

(a)  Effect  of  Subcooling 

Most  of  the  mechanisms  proposed  have  a  heat  flux  depend- 
ing strongly  on  the  subcooling.  This  is  not  the  case  in 
nucleate  boiling  as  revealed  in  many  experiments.  On  the 
contrary,  the  triple  interface  evaporating  mechanism  shows 
that  heat  flux  in  nucleate  boiling  depends  on  the  superheat 

At    =  f-t    only.  From  Eq.  (20),  although  bubble  pop- 
sat    w  sat 

ulation  N  and  diameter  d  also  play  important  roles,  the  bub- 
ble population  can  be  considered  as  constant  for  a  given  wall 
temperature  because  it  depends  only  on  the  temperature  of  the 
liquid  that  is  in  contact  with  the  wall.  However,  as  the  bub- 
ble diameter  decreases  with  the  increase  of  subcooling,  the 
heat  flux  contributed  by  nucleate  boiling  will  be  expected  to 
be  less  in  a  high  subcooled  liquid.  This  decrease  combines 
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with  the  increase  of  the  part  contributed  by  stirred  convec- 
tion will  give  an  observed  heat  flux  highly  insensitive  to 
the  degree  of  subcooling. 

For  a  given  pressure,  bubble  diameter  d^  is  smaller  in 
a  higher  subcooled  condition,  and  a  smaller  d^  will  make  a 
greater  At£at  since  more  bubbles  are  required  to  reach  the 
burnout  point  where  maximum  heat  flux  is  concerned.  Taking 
these  facts  into  consideration,  Eq.  (28)  shows  that  subcool- 
ing should  increase  (<lA)max» 

(b)  Effect  of  velocity 

In  the  region  of  strong  nucleate  boiling,  the  influence 
of  velocity  is  small.  This  can  be  explained  in  the  same  way 
as  in  the  preceding.  Since  the  convective  velocity  is  not 
able  to  change  the  temperature  of  the  liquid  in  contact  with 
the  wall,  bubble  population  N  will  remain  unaffected.  Fur- 
thermore, because  of  the  unstable  force  resulting  from  the 
velocity,  bubble  diameter  d0  is  going  to  be  decreased  as 
velocity  increases.  Then  the  heat  flux  reduction  due  to  this 
decrease  of  bubble  diameter  is  compensated  for  by  the  increase 
due  to  forced  convection.  This  again  will  make  the  total  heat 
flux  highly  independent  of  velocity.  Also  the  maximum  heat 
flux  can  be  increased  by  increasing  velocity  to  a  certain 
extent. 
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(c)  Effect  of  the  thermal  conductivity  of  heat  transfer 
surface  material. 

In  many  other  mechanisms,  heat  is  assumed  to  flow  through 
the  whole  part  of  the  heating  surface  that  is  not  covered  by 
bubbles  and  the  thermal  conductivity  can  be  generally  neglected 
in  comparison  with  the  very  low  conductivity  of  liquid.   But 
in  this  new  mechanism,  heat  flows  through  a  very  narrow  area 
near  the  triple  interface.   Thermal  conductivity  is  essential, 
as  can  be  judged  from  the  rapid  cooling  effect  observed  by 

Moore  and  Mesler  (24). 

Because  of  the  cooling  effect  of  triple  interface  evapor- 
ation, the  actual  (Atsat)act  (locally  near  the  bubble)  is 
different  from  the  observed  Atsat  (average  of  the  whole  sur- 
face). The  difference  Atsat  -  (^Sat>actual  is  Rurally 
dependent  on  thermal  conductivity.  Higher  conductivity  has 
a  lower  cooling  effect,  thus  a  small  difference  between  actual 
Atsat  and  observed  Atsat.  From  this  it  can  be  concluded  that 
for  a  given  observed  Atsat,  a  higher  thermal  conductivity 
will  give  higher  heat  flux. 

Since  the  effect  of  conductivity  only  shifts  an  observed 
temperature  scale  to  an  actual  scale,  the  maximum  heat  flux 
is  not  affected.   Of  course,  for  a  low  conductive  heater, 
more  observed  Atsat  is  required  to  attain  the  same  amount  of 
maximum  heat  flux. 
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(d)  Effect  of  surface  tension 

For  small  values  of  surface  tension,  the  bubble  is  ex- 
cited at  lessAtsat.   That  is  for  a  given  Atsat>  the  lower 
the  surface  tension  the  higher  the  heat  flux.  But  as  the  low 
surface  tension  liquid  reaches  its  burnout  point  in  a  smaller 

4t'  ...  a  lower  peak  heat  flux  would  result. 
sat 

(e)  Effect  of  pressure  and  temperature 

Although  the  effect  of  pressure  and  temperature  seems  to 
be  two  different  things,  in  the  triple  interface  evaporation 
mechanism  it  can  be  considered  as  only  one  factor.  This  is 
because  only  the  saturation  temperature  which  is  determined 
by  the  pressure  is  of  concern. 

The  effect  of  pressure  and  temperature  concerns  surface 
tension,  latent  heat  and  temperature.  For  a  given  liquid,  of 
course,  there  are  definite  relations  between  surface  tension 
and  temperature  and  between  latent  heat  and  temperature.  But 
these  are  only  in  experimental  forms.  Equation  (2S)  can  be 
further  simplified  by  substituting  some  typical  relations 
between  d^,  At'  t  and  some  known  parameters  such  as  surface 
tension,  specific  volume,  gravitational  acceleration,  etc. 
For  preliminary  treatment  of  the  new  mechanism,  however, 
Eq.  (2#)  will  remain  as  given  and  a  qualitative  discussion 
of  the  influence  of  pressure  and  temperature  will  be  given. 

From  Eq.  (20),  the  temperature  term 
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-T-   exp  ("  wr—] 

RTsat 

varies  exponentially  with  absolute  temperature.   Figure  22 

(A)  shows  the  case  of  water  whose  heat  transfer  coefficient 

at  1000  R  (962.5  psia)  is  known  to  be  200  times  greater  than 

the  value  for  the  same  At    at  atmospheric  pressure  if  the 

sat 

effect  of  N  and  dQ  is  neglected  temporarily.  Since  the  sur- 
face tension  decreases  as  temperature  increases,  the  bubble 
population  increases  while  the  detaching  bubble  diameter 
decreases  with  the  increase  of  temperature.  Experimentally 
the  variation  of  population  is  much  more  than  the  variation 
in  diameter.  Thus  if  the  effect  of  N  and  dQ  is  considered, 
a  still  higher  coefficient  will  result  in  the  high  temper- 
ature range. 

For  maximum  heat  flux,  since  the  surface  tension 
diminishes  gradually  to  zero  at  the  critical  state,  the 

At'   required  to  reach  the  peak  flux  also  diminishes  due 
sat 

to  the  ease  of  vaporization.  Qualitatively  the  product  of 
>/n"'  At'  .  can  be  expressed  as  in  Fig.  22  (B).  From  Eq.  (27) 
and  Figs.  22  (A)  and  (B),  it  can  be  seen  that  the  maximum 
heat  flux  in  nucleate  boiling  first  increases  with  the  in- 
crease of  pressure  until  a  certain  value  is  reached.   It 
then  drops  gradually  to  zero  at  the  critical  pressure. 

As  a  summary  of  the  above  discussion,  the  results  of 
these  conclusions  are  plotted  in  Fig.  23.  From  this  figure 
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it  is  clear  that  all  the  factors  that  affect  nucleate  boiling 
heat  transfer  are  explained  very  well  by  the  new  triple  inter- 
face  evaporation  mechanism. 


3-5.  Conclusion 


Although  many  excellent  works  show  that  bubble  stirring 
action  dominates  the  heat  transfer  rate  at  low  nucleate  boil- 
ing ranges,  there  are  facts  which  reveal  that  at  a  high 
nucleate  boiling  region,  heat  flux  is  mostly  contributed  by 
latent  heat  transport.  The  large  amount  of  heat  that  is 
transported  passes  through  the  bubble  via  a  very  small  area 
near  the  solid-liquid- vapor  triple  interface  at  the  bubble 
base.  The  heat  flux  is  entirely  determined  by  the  super- 

'  heat ,'  At     +   and  the  total  bubble  circumferences  that  are 
'   saw 

in  contact  with  the  heating  surface.  Prediction  of  heat 
flux  under  any  circumstance  will  be  possible  if  the  relation 
between  bubble  population,  bubble  diameter  and  superheat  is 

clear. 

From  an  available  relation  observed  by  Gaertner  and 
Westwater  (33),  the  triple  interface  evaporation  mechanism 
shows  that  the  heat  flux  is  proportional  to  N0*4   in  the 
high  flux  range  and  varies  directly  with  N  in  the  low  flux 
range.   Both  results  agree  with  the  experiments  of  Gaertner, 
Westwater,  and  Jakob. 

As  the  most  unexplainable  factors  in  most  proposed 
mechanisms,  subcooling  and  velocity  are  shown  in  the  new 


55 


mechanism  to  give  no  more  heat  flux  at  a  given  superheat 
Atsat,  while  they  do  increase  the  peak  flux  to  a  certain  ex- 
tent. As  a  preliminary  treatment,  a  single  equation  which 
successfully  explains  all  the  factors,  subcooling,  velocity, 
thermal  conductivity,  surface  tension  and  pressure  that  affect 
the  necleate  boiling  heat  transfer  has  been  derived  by  the 
author  of  this  report. 
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NOMENCLATURE 

A        ■  area  of  heat  transmission,  ftc 

a        =  evaporation  constant  defined  in  Eq.  (13), 

lb/hr  ft 

b        -  constant  defined  in  Eq.  (17),  ft2/hr 

d        =  instantaneous  bubble  diameter,  ft 

d.       =  instantaneous  bubble  contact  diameter,  ft 
c 

d        =  detaching  bubble  diameter,  ft 
o 

d»       =  detaching  bubble  diameter  at  critical  heat  flux, 
o 

ft 

d*       -  mean  effective  bubble  diameter  defined  in  Eq. 

(19),  ft 

=  thermal  conductivity,  BTU/hr  ft  deg  F 

=  Boltzmann  gas  constant,  BTU/molecule  deg  F 

=  latent  heat  of  vaporization,  BTU/lb 

■  mass  rate  of  evaporation,  lb/hr  ft2 

-  bubble  population,  Bubbles/ft2 
=  bubble  population  at  peak  heat  flux,  bubbles/ft 
=  pressure,  lb/ft 
=  vapor  phase  pressure,  lb/ft 
=  heat  rate,  BTU/hr 

-  gas  constant,  BTU/lb  deg  F 

-  bubble  radius,  ft 
=  bubble  radius,  ft 
«=  absolute  temperature,  deg  R 
=  absolute  saturation  temperature,  deg  R 


K 

k 

L 

m 

N 

Nf 

P 

pv 

q 

R 

R 

r 

T 

T„„„ 
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T 
w 

B 

absolute  wall  temperature,  deg  R 

Too 

- 

absolute  bulk  temperature,  deg  R 

V 

= 

volume ,  f t-* 

P 

= 

contact  angle,  deg 

I 

= 

effective  thickness  defined  in  Eq. 

(12),  ft 

6 

= 

surface  tension,  lb/ft 

? 

= 

density,  lb/ft3 

fr 

= 

0 

density  of  vapor,  lb/ft 

?L 

- 

density  of  liquid,  lb/ft ^ 

At 

= 

temperature  difference,  deg  F 

Atsat 

- 

VW  deg  F 

At'  ■  degree  of  superheat  at  peak  flux,  deg  F 
sat 

6  =  time,  hr 

"Jf  =  constant  defined  in  Eq.  (22) 

NNU  =  bubble  Nusselt  number,  dimensionless 

N  ■  Reynold  number,  dimensionless 
RE 

N  =  Prandtl  number,  dimensionless 

(q/A)LH  =  heat  flux  due  to  latent  heat  transport,  BTU/hr 

ft2 

(q/A)TrT  =  total  heat  flux,  BTU/hr  ft 
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APPENDIX  1.   CONVERGENCE  OF  EQUATION  (10). 
The  series, 


x»o 


becomes  infinite  when  y  =  0.  Since  at  y  -  0, 


r.O 


St  !     _.  4(tar-tSat)  y|_  ^ 

dX  lx=»  y=o      &      © 


But  for  any  y  >0  (y  *   0),  the  series  can  be  written, 

<a-y)C2ntQTT  ttL-y)f2nt-0TT 


•    e"      -  e 


4CU-W)  ^  '*. 


e-2ir 


<   1-e2"  V  e 


65 


.£5 

Since   Z  O**™       is  a  geometric  series  with  common  ratto  e 


2S  -tznn         .     ..„_._ .  „  ..j.,.1.  _™_  v.^4-,-0 


0 


for  any  y  =  5  ,   §  >0 


Moreover , 


p  <a  =  ! <•  i 

^  gny  ^  -1 


6       <   6 


00   -I-2J1TI 


e  a     converges 

•'  o 

from  the  Weierstrass  Comparation  Test,  it  is  known  that  Eq. 
(10)  converges  uniformly  in  the  interval  Z^y^CL,    %>0, 
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APPENDIX  2.      NUMERICAL  VALUATION  OF  Q5_aFROM  EQ.    (11). 
Equation    (11)   is  written, 
^-IX&jf^f,^  CM) 

Let    a-S=ma,    i.e.   m=l--|-   ,        where       0  6W1-<1,    Then, 
[l-Cosh^-^^'^JCosechCa^QTI 

=  [i-Cosh  wC2n+0ll]CosechC2n4-OTr 


.mCzm-Ott  ,    ^-rwftM+OTC 


mean*  0*1       e-wczn+uii 


_   r     __£ —  c    t    c? 1 

—  L    '         eC2rn-i)tt  _    p-(2n-H)Tt  j 

r      fi f2n+i)T[_       -f  Can^-OTI    j2 


(A) 


For  m  *0.5,  the  terms  e'^'^and  e"CWM*,)Tt  can  be  neg- 
lected for  all  n  except  n=0.  This  can  be  justified  by  cal- 
culating the  maximum  error  that  resulted  from  this  neglect, 
when  n=l,  m=0.5  as, 

[e°-75it_  q°-75K]  el,5rt 

e3TC    _     £-3IC  e3Tf 

6  =  [  e0,75K—  e~°-75n3a 

*   0.00934- 0-00905      any 
fc"      0.00334         C5'U* 

For  higher  n  and  m,  the  error  will  be  greatly  reduced;  thus 

this  action  is  reasonable.  Equation  (A)  becomes, 
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I- 


where  n  =  l,  2,  3,...,   m^O-5  Substituting  the  above  result 
into  Eq.  (11),  results  in 


Qs-a  = 


4K(L,-t«±)  {(e^-e^f 


IV  L  eir_  e-Tt  ,£,  (2M+I) 


fm-OC8n+i)3C 


Since 


rzm-o 


£0  (en+i)  x  S^M-X'' 


putting     eCm")Tr=x, 


gives 


The  value  of  TC  ®*"^; — r    for  various  m  is  tabulated  as 

4t<Ctw-tsot) 

follows: 


Yr\ 


o.o 

0.5 

0.7 

0.9 

0.95 

0.99 


s 
a 

KQs-cl 

4Kft«,-t»0 

1.0 

0.000 

0.5 

0.135 

0.3 

0.332 

0.1 

0.345 

0.05 

1.172 

0.01 

1.844 
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APPENDIX  3.   EFFECTIVE  THICKNESS  OF  WATER. 

In  Section  3-2  it  is  shown  that  for  water  at  peak  heat 
flux,  the  total  heat  transfer  area  is  20$  of  the  heating  sur- 
face  area.  Thus  for  a  semi-spherical  bubble,  Eq.  (26)  gives 
a  heat  transfer  area  per  unit  heating  surface  area  as 

TURN'S  =  TT*d.'x  1.15  Cf-jj)-  8-15  S-^r 

Since  this  area  should  equal  0.2,  that  is 

8.15  S  =  0.2  d'o 

do' 


or 


S 


40.75 


which  shows  the  effective  thickness  is  really  very  small, 
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ABSTRACT 


Since  I960,  the  latent  heat  transport  in  nucleate  boiling 
heat  transfer  has  received  considerable  study  and  its  contri- 
bution has  been  revised.  Recently  it  has  been  observed  by 
several  investigators  that  latent  heat  transport  is  signif- 
icant at  all  stages  of  nucleate  boiling.  The  purpose  of 
this  report  is  to  further  verify  the  contribution  of  latent 
heat  transport  and  to  give  a  theoretical  analysis  of  this 
mechanism. 

The  data  of  Westwater  and  Santangelo  show  that  the 
latent  heat  transport  contributes  about  50%  of  the  total 
heat  flux  observed  in  boiling  saturated  liquid  at  moderate 
heat  flux  range.  Careful  examination  reveals  that  the  heat 
flow  into  the  bubble  from  heating  surface  passes  through  a 
very  narrow  area  near  the  solid-liquid-vapor  triple  inter- 
face. A  triple  interface  evaporation  mechanism  is  thus 
proposed  on  the  basis  of  this  investigation.  The  equation 
of  bubble  growth  rate  derived  by  this  mechanism  has  the  form 


R  =  Const  J  6  tan  -|- 

This  agrees  with  the  recent  observation  made  by  Johnson,  Jr. 
and  coworkers.  The  correlation  of  heat  flux  contributed  by 
latent  heat  transport  as  derived  by  the  author  of  this  report 
is 

q/A  =  Const  N  dQ  L2   exp  (  L )  Atgat 

RTsat      RTsat 


This  equation  has  been  used  successfully  to  explain  the 
mechanism  by  which  the.  factors  pressure,  velocity,  degrees 
of  subcooling,  surface  tension,  and  thermal  conductivity  of 
the  heating  surface  affect  nucleate  boiling. 


